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Abstract 

Employing a simple calculation method obtained by M.-H. Kato, we 
calculate the three loop renormalization group in the su(2) coset confor- 
mal field theory with a slightly relevant perturbation and the su(2) Wess- 
Zumino-Witten model with a particular invariant marginal perturbation. 
Zamolodchikov's c-theorem, exact data of the perturbation operator and 
a known exact form of the operator product coefficient enable us to cal- 
culate the beta function, the gamma function and the c-function to three 
loop order. This result gives the logarithmic finite size correction to the 
ground state energy and the low temperature behavior of the specific heat 
in the Heisenberg antiferromagnetic chain with high accuracy. We describe 
the consistency with results obtained by several authors on the basis of its 
exact solvability. We discuss an experiment of the specific heat and the 
suceptibility recently observed. 
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1 Introduction 



There have been extensive studies in conformal field theory (CFT) with a relevant 
perturbation which reveal universal nature of critical phenomena in low dimensional 
statistical physics. Landau- Ginzburg theory gives a good qualitative picture for the 
renormalization group (RG) flow including multi-critical points 0. In this theory, 
we can specify the driving relevant operator of the CFT at the ultraviolet (UV) 
fixed point and the irrelevant operator of another CFT at the infrared (IR) fixed 
point. We can check the expected flow practically in the e-expansion by calculating 
the deviation of the central charge which obeys Zamolodchikov's c-theorem 0. 

In this paper, we utilize Zamolodchikov's c-theorem and the exact form of the 
operator product expansion (OPE) coefficient to calculate the beta function to three 
loop order in the ^-expansion for su(2) coset models with the slightly relevant per- 
turbation with the dimension 2 — e. M. -H. Kato argued this application of the 
c-theorem to calculate the beta function up to two loop order for the slightly per- 
turbed general coset models in all A-D-E classes [[|. By taking the limit e — > 0, the 
coset CFT becomes a certain Wess-Zumino-Witten (WZW) model with a particular 
invariant marginal perturbation. In this limit, he also calculated the logarithmic 
finite size correction to the ground state energy and the logarithmic temperature 
dependence of the specific heat of quantum spin chains to two loop order. The ob- 
tained fitting function agrees with the groud state energy of an su{2) spin chain by 
a numerical Bethe ansatz [Q. Here, we calculate higher order logarithmic correction 
to the ground state energy and discuss the consistency with results obtained by sev- 
eral authors in different methods. This paper is organized as follows. In section 2, 
we review M.-H. Kato's method to calculate two loop beta function in general coset 
models with generalizing some equations for higher order calculation. In section 3, 
we confine our discussions to su{2) case. In su(2) coset model, the known exact 
form of the OPE coefficient of the driving operator enables us to calculate the three 
loop beta function. In section 4, we discuss the logarithmic finite size correction 
to low lying energy levels. We discuss consistency with a numerical Bethe ansatz 
and some other methods. In section 5, we calculate a logarithmic temperature de- 
pendence of the specific heat per unit length in a su{2) spin chain. We discuss the 
recent experiment [H from the view point of the higher order renormalization group 
calculation. 



2 Two loop calculation in general coset models 

Here, we review M.-H. Kato's method for two loop beta function Q with preparing 
some extended forms of the equations for three loop calculation in the next sec- 
tion. In this method, we consider perturbed coset CFT A4(k, I; G) = Gk® Gi/Gk+i 
by a slightly relevant operator $uv with self-closing algebra which corresponds to 
the (1,3) operator in the Virasoro minimal model clS cl TQ gularized theory for the 
marginally perturbed WZW model. In this section, we assume those following two 
facts obtained by Ahn, Bernard and Leclair [[F]. They showed that the slightly 
relevant operator <3>uv drives the CFT at UV fixed point to the IR one with the 



irrelevant perturbation of the operator as in su(2) case [[I], 



l)+gf ^w(z, z) - l-k)+g'J — r<M^), (IR limit) (1) 

where H(k,l) is a critical hamiltonian in Ai(k,l; G). They showed also that this 
deformed CFT A4(k, I; G) by the operator <3>uv becomes the WZW model with 
an invariant marginal perturbation, that is 

G k © Gi/Gk+i — > Gfc, (/ — > oo), 

and 

$ uv (z, f) - f) = --4= £ ■'"CO^®. (' - oo), (2) 

in a limit I — > oo, where Z) is the dimension of the Lie algebra G. The existence of 
RG flow between these fixed points is guaranteed by the self-closing algebra of the 
operator $uv- m the CFT Ai(k, I; G), the central charge is 

c(k,l;G) = c(G k ) + c(Gj) - c(G k+l ) (3) 
kD ( x h c (k + h c ) \ 



k + h c \ (I + h c )(k + I + h c ) J 
and the conformal dimensions of the operators $uv and $ir are 



Am = 1 + 



I — k + h c 



where r is the rank and h c is the dual coxeter number of the Lie algebra G. The 
dimension of the Lie algebra is given by D = r(l + h c ) in simply- laced algebra. For 
example in G = su(n), these are given by r = n — 1, h c = n and D = n 2 — 1. Here, 
we define the parameter 

and discuss the RG flow from the UV theory A4(k, I; G) to the IR theory Ai(k, I — 
k; G) in the e-expansion. The limit e \ as I /* oo gives 

Auv = 1 - | - I- 

This relation is necessary for the identification of the operator (pD, which should 
be shown by the correspondence between these OPE coefficients with their closing 
algebra, as well as their conformal dimensions. The OPE relations of the operator 
$uv and the marginal operator are 

$uvOM)$uv(0,0) ~ Ji|L$ uv ( ,0) (6) 

$(M)$(0,0) ~ ,^$(0,0). (7) 



with b(s) — > bo as e — > 0. To obtain the OPE coefficient b(e) in the e-expansion, we 
assume Zamolodchikov's c-theorem |2|. The beta function of the deformed CFT by 
the operator $uv with the running coupling constant g is 

m = - Eg + b M 9 ' + d M g * + eM g * + ... : (8) 

and the gamma function of the operator $uv is given by 

7( 9 ) = 2+^. (9) 
The coefficients in the beta function is expanded in e 

oo 

b(e) = (10) 

n=0 
oo 

d(e) = X>n£ n , (11) 

n=0 

oo 

e(e) = J2 e ne n , (12) 

n=0 

(13) 

The beta function @ has a trivial fixed point g> = and another non-trivial one 
g = g* ^ 0. The trivial fixed point g = is the UV CFT Ai(k, I; G) and the other 
one corresponds to the IR CFT A4(k,l — k;G). The non-trivial fixed point is 
expanded in a series 



9* = E 9n£ n , (14) 



n=l 



whose coefficients g n are written in terms of those of the beta function (O) 

= r« (15) 

2 



6§ 

f3 = i (oo°i - °o°2 + 6o Mo + 8dg - 26gdi - 46 e ) 



6g 

Zamolodchikov's c-function is defined by 



(16) 



6 [9 

c{g) = cuv + - / P(x)dx, (17) 

where cuv is the central charge of the CFT A4(k, I; G). Zamolodchikov's c-theorem 
gives us the following constraints 

cm = c{jg m ), (18) 
A IR = 7 (^)/2 = l+/3'(^)/2, (19) 



where cir is the central charge of the IR CFT Ai(k, I — k; G). The right hand side 
in these constraints can be expanded in e 

e 3 {2b! 3d \ , ( 2b 2 3di - 36? -60b lC l + 24e 12d \ = 
c m -c„ v = -jj +^ + e + [^ + -jg— + jjj -jj-j e +■ 

(20) 

and 

Air = 1 + £ - + % 2 + i (-26 Mo - 4dg + &J& + 46 e ) e 3 + • • • . (21) 

The knowledge of the flow Ai(k,l;G) — > A4(k,l — k;G) by the slightly relevant 
perturbation gives further constraints 

- - (22) 

The £-expansion of these constraints (^) and (El) can be done as follows: 

CIR - CUV = __^ 1 + _ e + _ e > + ...j | (24) 

a, r = i + £ + A e2 + J £ 3 + .... (25) 

These constraints (|20|) , ([2l"l), (PH) and (|25| ) can be solved with respect to b\ and c? 
in terms of b 

where &o is OPE coefficient of the marginal operator in the WZW model 

2h r 



k^D' 



Therefore, we obtain the beta function of the marginally perturbed Gk WZW model 
up to two loop order 



without explicit calculation. 



3 Three loop beta function 

In the previous section, we obtained the two loop beta function with the help of two 
assumptions ([!]) and (|j). In this section, we limit our discussion to G — su{2) case 
and we check these assumption within the framework of the ^-expansion. In this 
case, we have D = 3, h c = 2 and r = 1. We have the slightly relevant operator $uv 



with the conformal dimension Auv — 1 — 2/(fc + / + 2). The exact form of its OPE 
coefficient b(e) in the CFT M.(k, I; su(2)) is given in || 

h(*\ - 8 (2-{k + l)e) 2 
{) ~ V3k(2-e)(4-(k + 2)e) [Z ° 



x 



r(i-g) / r(i- g /4) \ 3/2 /rxi + 3£/4) \ 1/2 / r(i + g /2) \ 2 

r(l + e) W+e/4); {v{l-3e/A)J \T{l-e/2)J 



This gives us the data of b , b\, b 2 , ■ ■ ■ and then the method explained in the pre- 
vious section enables us to calculate the beta function up to third order. The OPE 



coefficient ( j28|) is expanded in 



4 / 3k k (fc-6 ) , \ 

b(e) = —=- 1 e + — -e 2 + ■■■). 28 

1 ; V3k V 4 16 J 1 ; 

This formula shows 

b{£) ^ b ° = ^kk ( ^ 0) ' 

which justifies the identification (^|) we assumed. The one loop beta function is 
given by /3(g) = bog 2 / 2 which enables us to calculate the central charge cir and 
the conformal dimension Air by eq. (|20|) and eq.(^). These results show that the 
expected flow (Hj) is checked in the ^-expansion around the UV fixed point at one 
loop level. Therefore, we can use the constrains to solve the higher order coefficients 
in the e-expansion. Since the constraints (|2"0), (|2T|), (23) and ( p5|) can be solved 



with respect to do, d\, and eo in terms of bo, &i, and 62, we obtain them explicitly 

4 10(fc-2) 10(fc-2) 

«o = — r, d\ = : , e = 7= , 29 

3fc' 1 9k 9V3k 2 V ; 

where the rank and the dual coxeter number of su(2) are given by r = 1 and h c = 2. 
The limit e — > gives the beta function to three loop order for the marginally 
perturbed level k su{2) WZ W model 

m = b j9 2 + ^9 3 + e j9 4 + 0(g% (30) 



4 The c-function and the logarithmic correction 

Here, we consider the perturbed WZW model in the asymptotically non-free re- 
gion g > 0. The marginally perturbed level k su{2) WZW model in asymptotically 
non-free region describes S = k/2 antiferromagnetic Faddeev-Takhtajan model. Es- 
pecially for k — 1, this corresponds to S — 1/2 Heisenberg antiferromagnetic chain 
model. We calculate the c-function by integrating the beta function and discuss the 
logarithmic finite size correction to the ground state energy of the antiferromagnetic 
chain. First, we integrate the differential equation 



(31) 



for the running coupling constant g(l) 



I{g{l)) -I(g(a)) = log l/a, (32) 

where a is the UV cut off length scale of the field theory. Here, we define the integral 
function I{x) by 

f dx 2 2d , 2 / dl \ „ . 9x . , 

Let L be the size of the chain which sufficiently large compared to the UV cut off a. 
In this case, we have g(L) << g(a) because of the asymptotic non-freedom. Here, 
we introduce the RG invariant length scale L = aexp(—I(g(a))). The running 
coupling constant g(L) is determined by 

I(g(L)) = log L/L , 

and solve this equation by expanding l/I(g(L)) in g(L) and s = l/log(L/L ) 

,(L) = ^(l + .| bg . + ^(^.-*(*±M) + ...). (34) 
Substituting this running coupling constant into the c-function (0), we obtain 



c(g(L)) = c + 3 + 4 (35) 

log L/Lq log £/£ 



A 3 + A' 3 log log L/L + A 3 '(log log L/L 



>2 



+ J " -^.^^u ^^^^-w + O(log- b L/L ), 
log L/L 

where the central charge cq of level k su(2) WZW model and other coefficients are 

Co = why (36) 

A l = -k\ 

8 ' 

q q 

64 ' 2 16 

q q q 

16 3 32 3 16 

This formula has only one fitting parameter L Q . The logarithmic finite size correction 
to the ground state energy E (L) in the Faddeev-Takhtajan chain with finite length 
L is calculated by 

E (L) = eoo L - —c(g(L)) (37) 
6L 

with the formula (pB]), where is energy density in the infinite length limit and v 
is the spin wave velocity. This formula ([36]) fits well the data of the ground state 
energy in a numerical Bethe ansatz obtained by Nomura ||] in k = 1 case. He shows 
the row data of low lying energy levels in the Heisenberg antiferromagnetic chain 
with L = 256 - 16384 lattice sites with unit lattice spacing. To check the consistency 



with all data including the excited energy, we have to consider the gamma functions 
of two primary fields in WZW model. In the level one su{2) WZW model, we 
have four primary fields which correspond a singlet and a triplet excitations with 
conformal dimension (1/4, 1/4). We write their gamma function 

i\g) = 7S + 7^ + ^(M 2 + 7^(M 3 + 0(^ 4 ), (38) 
l S (g) = 7o + 7?M + 72(M 2 + 73(M 3 + 0(^ 4 ). 

The level one su(2) WZW model gives 

7o = 7o = 1/2, 

and one loop RG calculation gives 

7| = -1/2, 7i S = 3/2, b = ^=, 

however other coefficients have never been calculated. Note that the first excited 
state has three fold degeneracy because of the su{2) symmetry. First, we fit the data 
at two loop level. In this order, the invariant scale is given by L Q 1 = 1.74087 for 
fitting the data of the ground state energy in J| by the formula (^). Unfortunately 
in the present stage, we have only the first order gamma function of the primary 
fields. Though we cannot check his data of excited energy in second order explicitly, 
we can fix the coefficients 72 and 7I by fitting the data as follows: 



1 1. 
1 

2 2 

1 3, 



<Loo P (g) = 1-00012 (-- -b g) - 0.0619577(M 2 , (39) 



l2Loo P (g) = 0.999499 (J + |MJ + L00147(M 2 - 

Next, we discuss three loop order. The invariant length scale becomes L 1 = 2.14538 
to fit the ground state energy in this order. The gamma functions for fitting the 
excited energy become 

ikoop(g) = 0.999868 Q - - 0.00772577(M 2 - 0.284163(M 3 , (40) 

lkoo P (g) = 0.999932 (i + ^&os) + 0.863347(^) 2 + 0.64291(^) 3 . 

It is important that both coefficients of the first order are quite close to unity as in 
two loop order. This fact shows a good agreement between our three loop calculation 
and the numerical Bethe ansatz. On the other hand, other coefficients change from 
those in the two loop order. 

In this fitting, we compare our calculation of some typical values with those by 
Nomura in the numerical Bethe ansatz 

c Be the{g{256)) = 1.00103233, c Bette (g(16384)) = 1.000239164, (41) 
C2Loo P (g(25Q)) = 1.00101926, c 2ioop (#( 16384)) = 1.000241815, 
C3Loop(g(256)) = 1.00101115, c 3 Loop(0(16384)) = 1.000243613, 



7^0/(256)) 

liLooM™)) 



0.464830, 
0.464835, 
0.464832, 



7^(^(16384)) 
7^0/(16384)) 
7^(^(16384)) 



0.478301, 
0.478296, 



0.478297 



(42) 



7WS(256)) 
7 2 s l oop (^(256)) 
7 3 s Loop (^(256)) 



0.620490, 
0.620513, 
0.620493, 



7^(^(16384)) 

7 2 W2(16384)) 
7 3 \ oop (<?(16384)) 



0.570916, 
0.570910, 



0.570912 



(43) 



The two loop calculation is slightly better than the three loop one with respect to the 
c-function. It is not less trivial that this three loop gamma functions fit the excited 
energy levels very well, since this fitting has two more parameters 73,73. These 
results of the RG calculation seem consistent with the numerical Bethe ansatz. 

Here, we discuss consistency of our results with those obtained by several authors 
0, |9], [K| [LlJ]. First, we comment on Nomura's fitting of his numerical data based 
on a two loop RG calculation M, which gives A\ = 0.365162. He take into account 
only first order of corrections to the c-function, although he uses the running cou- 
pling constant obtained from the two loop beta function. This is consistent way 
to take into account log log L/Lq log -4 L/Lq term in the c-function (|3~6|). However, 
we cannot neglect log~ 4 L/Lo term to fit the numerical data, since log _4 L/Lo term 
behaves as large as log log L/Lq log -4 L/Lq in numerical data of the model with not 
so large degrees of freedom (256-16384 sites). To obtain more accurate fitting, we 
should take into account g 4 order in c-function and g 2 order in the gamma function 
of the primary fields. He fixes the parameter Lq 1 = 0.56532 in such a way that 
the one loop correction to the triplet excitation energy with the two loop running 
coupling constant fits the numerical data of the triplet excitation. It is possible to 
fit only one observable by adjusting the parameter Lq with neglecting the second 
order correction in it, since the second order logarithmic correction can be absorbed 
into the first order one by adjusting the parameter Lq. In this procedure however, 
the second order logarithmic correction cannot be absorbed into first order one to fit 
observables more than one without the second order correction. Despite his recogni- 
tion of the g 2 order correction to the singlet excitation energy and L~ 2 correction to 
the ground state energy by the irrelevant operator L_ 2 X_ 2 1, his obtained c-function 

is c(g(L)) = 1 + 0.36516 + ^jfr-, which gives not so good agreement with 

the leading coefficient A\ = 0.375. If we take care of log _4 L/L term from the g A 
term in the c-function and log -2 L/Lq term from the g 2 term in the gamma function 
7*(g), the numerical data should be fitted with sufficient accuracy even in two loop 
order. Actually, Karbach and Mutter obtain A\ = 0.375 A' 2 = —0.73 and A 2 = 0.15 
with recognition of log _4 L/L term in the c-function with a careful treatment of 
fitting their data in numerical Bethe ansatz. Their obtained A% agrees with ours, 
and also Kliimper shows this result A x = 0.375 by calculating the specific heat in 
his new method of solving a set of non-linear equations based on a quantum transfer 
matrix for the thermodymamic quantities partially using a numerical calculation. 
We will see the relation between the specific heat and the ground state energy in 
the next section. Next, we discuss the two loop order coefficients. Our results 
A' 2 = -9/16 = -0.5625 and A 2 = 9/64 + 9/8 logL ~ -0.483058, with the two loop 
fitting parameter Lq 1 = 1.74087 in (j36D does not agree with those by Karbach and 



Mutter. We expect our formula with the three loop correction ( |3"6| ) will be useful to 
fit their data, eventhough fitting the numerical data for the higher order coefficients 
seems extremly difficult. 

Finally, we mention the results obtained by Woynarovich and Eckel with Euler- 



MacLaurin expansion in analytic Bethe ansatz [] I ] . The coefficients 7q, 7*, 7q and 



7f in excited energy level and Co in the WZW model agree with those obtained 
by Euler-MacLaurin expansion in the Bethe ansatz. On the other hand, the first 
coefficient in the logarithmic correction in eq. (|37|) is calculated as A\ = 0.343347 ■ • • , 
in |TT] which does not agree with the result A\ = 3/8 = 0.375 obtained by the WZW 



model. This is well-known contradiction pointed out in [|| ||, [10|, 12| . Since our fitting 



of numerical Bethe ansatz, an appropriate treatment of the numerical Bethe ansatz 
and Kliimper's result [[Kj] support the WZW model, we believe A 1 = 3/8. 



5 Specific heat and suceptibility 

Next, we calculate the specific heat per unit length. The partition function of one 
dimensional quantum system is given by 

Z(L,M) = tre- MH/hv . (44) 

M = has the dimension of length, where v is a spin wave velocity, k-Q is the 
Boltzmann constant and T is the temperature. The partition function possesses a 
modular invar iance 

Z(L,M) = Z(M,L), (45) 

and the free energy per unit length f(L,M) = —^-\ogZ(L,M) possesses it as 
well. In the low temperature limit M — > 00, the free energy becomes 

/(L,oo) = E (L)/L (46) 
= - — c (g(L)). 

The modular invariance gives the temperature dependence of the free energy per 
unit length for the infinitely long chain as follows: 

/(oo,M) = f(M, 00) (47) 

The specific heat per unit length C{T) = -Tf£ is calculated in a low temperature 
expansion 



Shv log Tq/T log Tq/T 



B 3 + ^loglogT /T + ^(loglogTo/T) 2 6 

+ l^njf +0(log To/T) 



where To = -j^j^ is a parameter with the dimension of temperature. The coefficients 
are given by 

Co = 1+2 (49) 



B l = -k 2 



B 2 = —k 2 (k + 4), B'=- — k s , 
2 64 v ; ' 2 16 ' 

B 3 = -—k 2 (k-4), B'=-—k 3 (k + A), B" = —k 4 . 
This formula shows Ai = A' 2 = _B 2 and A3 = 5 3 ' given in eq.(^), and therefore 



our result is consistent with that obtained by Kliimper for the specific heat [[10 
The specific heat per unit length of any model in the universality class including 
S = 1/2 Heisenberg antiferromagnetic chain obeys this result with k = 1. For 
example, the specific heat per unit length of S = 1/2 antiferromagnetic chain with 
a small perturbation of a second neighbor interaction satisfies the same formula by 
adjusting To. This shift occurs through the change of the initial coupling constant 
g(a). As proposed in |J, this low temperature behavior of the specific heat per 
unit length (|49| ) in the S = 1/2 (k = 1) Heisenberg chain will have to be observed 
experimentally, as well as its suceptibility. Eggert, Affleck and Takahashi calculated 
its suceptibility by one loop renormalization group and the Bethe ansatz for the 
Heisenberg chain The result has been observed by N. Motoyama, H. Eisaki 

and S. Uchida, experimentally in an antiferromagnetic chain S^CuOs 0. Here we 
discuss the data fitting of both the suceptibility and the specific heat which are 
observed. The Heisenberg hamiltonian is 

H = J ^ Sj • Sj+i. 

i 

The spin wave velocity is given by 

v = nJ/2, 

which is determined from the Bethe ansatz. Eggert et al. gives the best fit To = 7. 7 J 
of the one loop renormalization group calculation for the suceptibility by the Bethe 
ansatz calculation of the Heisenberg chain [|13|]. They argued that the log -2 T /T 
term in the suceptibility can be eliminated by shifting T . In addition to the 
log -2 T /T term however, we have further logarithmic corrections from the two loop 
running coupling constant (|34D in the one loop corrected suceptibility 



(rp\ 1 1 (, loglogT /T\ n , _3 Tm 

x(T) = 2^ + 4 7 r,logTo/T ( 1 - 21ogT /T J + ° (1 ° g T ° /T) ' 

This formula should be better for fitting the suceptibility calculated by the Bethe 
ansatz. To fit the specific heat at the same time, however, the adjusted To for the 
suceptibility cannot be used for the specific heat with the same accuracy level as 
the suceptibility. This is because the same reson as in the fitting of numerical Bethe 
ansatz for the logarithmic finite size correction we argued in the previous section. 



The elimination of both the log -2 Tq/T term in the suceptibility and log -4 Tq/T term 
in the specific heat cannot be done by shifting the only one parameter at the same 
time. We have to calculate the coefficient of log~ 2 To/T term in the suceptibility to 
fit both specific heat and the suceptibility. To do so, we have to calculate the g 2 
order correction in the suceptibility explicitly. 

6 Conclusion 

In this paper, we have shown a simple calculation method for the higher order beta 
function of a perturbed CFT by utilizing the constraints given by Zamolodchikov's 
c-theorem. In a particular CFT model su(2)k © su(2)i/su(2)k+i, the known exact 
form of the OPE coefficient of the slightly relevant operator enables us to calculate 
the beta function up to three loop order in the e-expansion, where e = 4/(fc + 1 + 2). 
In the limit e —>■ as I — > oo, we have the level k WZW model with a marginal per- 
turbation. The obtained three loop beta function is useful to study a quantum spin 
chain model described in the marginally perturbed WZW model. The logarithmic 
finite size correction to the ground state energy in the su(2) quantum spin chain 
model has been calculated up to the log" 5 L order as well as the logarithmic temper- 
ature dependence of the specific heat. The obtained formula of finite size correction 
fits the data of a numerical Bethe ansatz quite well. The low temperature behavior 
of the specific heat obtained here should be observed experimentally. We provide a 
slightly better fitting function of the suceptibility than the one loop renormalization 
group calculation. Completely consistent treatment of both the suceptibility and 
the specific heat will be reported soon. 
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